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TECHNICAL MEMORANDUM 


1 . 0 Introduction 

Many substances are vented from various parts of 
Sky lab, and whenever this occurs a thrust results, which 
unless it is directed through the vehicle mass center, 
imparts a torque to the vehicle. The magnitude, direction, 
and duration of these torques determines the extent to 
which they can be nulled by the control moment gyroscope 
(CMG) attitude control system. This memorandum provides 
data on magnitudes, directions and time histories of 
significant Sky lab venting torques, and twenty-four hour 
profiles are included for use by those assessing CMG 
performance . 

2 . 0 Summary and Conclusions 

The Sky lab vents mentioned below give rise to 
propulsive effects of such magnitude that attitude control 
system performance may be impaired should one CMG fail. 

Waste Water Vent - Water accumulates as the CSM fuel 
cells are used, and if a full tank of water is dumped, 
an angular impulse of about 5000 ft- lb-sec results. 
Although the mission baseline calls for fuel cell 
shutdown after solar array deployment, the option to 
continue operation until reactant depletion is being 
kept open until launch. At a 1200 watt power level, 
a full tank of water would accumulate in four days, 
but by dumping the tank more often, for example daily, 
the total impulse could be distributed more evenly. 

K„ and 0 ^ Relief Vents - With the fuel cells shut down, 

reactants must be vented from the cryogenic tanks. 
Initial angular momentum accumulation for the nominal 
mission is about 900 ft- lb-sec/orbit for H 2 and half 

this amount for C> 2 . These values decrease as the 

tanks empty? thus, momentum/orbit will decrease with 
time, and the initial value will decrease if delays 
are experienced in attaining orbit. The C> 2 value 
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above (450 ft- lb-sec/orbit) arises from the oxygen 
expelled from the tanks in excess of that required 
for metabolic and cabin leakage makeup. If a temporary 
increase in cabin C> 2 partial pressure - about 1/4 psi - 

is allowed, 0 2 venting would be unnecessary. 

LOX Tank Vents - If a day's collection Of urine (13 lb) 
enters the tank in liquid form, about 1200 ft-lb-sec 
will be accumulated in the first orbit with smaller 
amounts in the five succeeding orbits. A trash dump 
containing four pounds of water results in about 800 
ft-lb-sec in the first orbit. The above values rise 
to 19,000 and 13,000 ft-lb-sec, respectively, if one 
of the two "nonpropulsive" vents should freeze shut. 

(A heater will be placed on one vent to avert this 
possibility. ) 

M092 Experiment - This vent lies under the OWS meteroid 
shield, and the flow will be diffused somewhat. With 
no diffusion, about 1300 ft-lb-sec accumulate each time 
the Lower Body Negative Pressure Experiment is performed. 
MSFC estimates that diffusion will reduce the magnitude 
to one third of the above value. 

That benefits can be obtained from scheduling and 
choice of options is illustrated by comparing the "typical" 
and "worst-case" day momentum profiles (Table IV) . Angular 
impulse imparted during most orbits of the typical day is 
about 200 ft-lb-sec, whereas on the worst-case day most 
orbits range from 40 to 900 ft-lb-sec. 

The data presented here are being used in a 
simulation of the Skylab CMG attitude control system to 
see if it performs satisfactorily with only two CMGs in 
operation.* Results from this work will be available shortly. 

It appears that- all typical day orbits and most worst- 
case day orbits can be accommodated. 

3 . 0 venting Sources Considered 

Table I, compiled basically from Reference 1, 
gives the significant vents, their Orbital Assembly co- 
ordinates, the direction of thrust, the effluent, and the 
maximum flow rate. Several comments concerning this table 
are in order. 
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THESE ARE CALCULATED RATHER THAN GIVEN VALUES. THE METHOD OF 
CALCULATION IS GIVEN ELSEWHERE IN THIS MEMORANDUM 
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The first three lines (waste water vent and the purge 
vents) represent one mode of operation in which all 
three fuel cells operate at a 1200 watt level until 
reactants are depleted (about 10 days) . The fourth 
and fifth lines (relief vents) represent an alternative 
mode in which fuel cells are shut down twenty-four 
hours after docking. In this case the remaining 
reactants are expelled from the tanks and vented. 

The choice of mode will be made before launch. 

The molecular sieve vents are nominally nonpropulsive ; 
the upper signs refer to one vent of the pair, and the 
lower signs refer to the other vent.* 

The urine dump and trash dump systems empty into the 
LOX tank, which then vents to space through a pair of 
so-called nonpropulsive vents. Again, the two sets 
of signs refer to the two opposing vents. 

4.0 Results for Individual Vents 


Results are summarized in Table II, to which the 
following comments apply. The direction cosines apply to 
either the torque or the angular momentum magnitudes given 
in the columns to the right. The torque given is the 
largest that will occur during the mission. Likewise, worst 
day max H/orbit is the largest angular impulse that will 
occur in one orbit of the mission. Typical day max H/orbit 
denotes the largest angular impulse that will occur in one 
orbit of a "typical" day. Note that there may be orbits in 
either the typical or worst days in which the impulse is 
much smaller than the max value stated. 

4.1 CSM Fuel Cells 


Two options are considered here. Option 2 is base- 
line, but option 1 will be kept open as an alternative until 
launch . 

Option 1 - Fuel Cells Operate 

Liquid water vents from the CM waste water vent 
under 25 psia of bladder pressure.** Thrust and flow rate 
may be obtained from Appendix A, formulas (A-6) and (A-9) , 
respectively, and Table II gives the resulting torque and 


* Data obtained from Ed Carmodv, McDonnell Douglas-East, 
November 2, 1970. 


**This figure obtained from Roger Tanner, MSC, June 26, 1970. 
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TABLE II 

RESULTS FOR INDIVIDUAL VENTS* 


VENT NAME 

DIRECTION COSINES 

MAX T 
FT-LB 

MAX H/ORBIT 
TYP. DAY** 

MAX H/ORBIT 
WORST DAY** 

X 

Y 

2 

CM WASTE WATER VENT 

.012 

-.993 

-.120 

3.5 

1130 

4530 

SM FUEL CELL H : PURGE VENT 

—.058 

-.010 

.998 

1.5 

NEGL 

.SGIBLE 

SM FUEL CELL 0 : PURGE VENT 

.053 

-.091 

-.994 

3.2 

NEGLIGIBLE 

SM H : RELIEF VENT 

-.058 

-.010 

.998 

0.17 

210 

930 

CM 0 : RELIEF VENT 

-.030 

-.866 

-.440 

.076 

0 

430 

AM CONDENSATE VENT 

-.138 

.882 

.450 

.042 

NEGL 

IGIBLE 

AM MOL. SIEVE VENTS 

.405 

-.749 

-.524 

.0072 

40 

40 

OWS URINE DUMP SYSTEM 

-.576 

-.797 

-.180 

0.94 

1160 

1160 

OWS TRASH DUMPS 

-.576 

-.797 

-.180 

0.70 

790 

790 

M092 LOWER BODY NEC. PRES. EXPT. 

-.051 ! 

-.776 

-.628 

1 A 

420 

1270 

M479 ZERO-G FLAM. EXPT. 

.103 ! 

.335 

-.937 

1.4 

0 

510 


• COMPUTED RELATIVE TO A MASS CENTER LOCATION OF (646.7, -2.1, -25.3) IN. IN ORBITAL 
ASSEMBLY COORDINATES 

** FT LB-SEC/ORBIT; 1 ORBIT = 5600 SEC 
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angular momentum accumulation. The waste water tank holds 
92 lb, and at the assumed power level, 20 lb/day of water 
are generated. We assume that typically water is dumped 
every day, and in the worst case, it is allowed to accumulate 
for four days, the respective dumping times being 320 sec 
and 1280 sec. 

Thrust for the fuel cell purge vents is obtained 
by the method of Reference 2. Due to the short purge intervals, 
the angular momenta associated with these vents is negligible 
compared to that of the waste water vent. 

Option 2 - Fuel Cells Shut Down 

Hydrogen is vented from the cryogenic storage tanks 
at a rate that depends on how full the tanks are (see Appendix B) . 
The peak flow rate, which occurs at fuel cell shutdown, is used 
for the worst-case day in Table II, and a 56-day-average flow 
rate is used for the typical day. 

Similarly, oxygen is vented from the storage tanks. 

Most of this 0 2 is used to replenish the cabin atmosphere, and 

the excess, if any, is assumed to be vented overboard. Flow 
rate calculations are covered in Appendix B; the worst-case 
day again occurs at fuel cell shutdown, but since overboard 
venting lasts for six days at most, the typical day flow rate 
is zero. 

4.2 AM Condensate Vent 


Thrust is calculated as shown in Appendix A. Onlv 
about five minutes is required to dump a three-day accumulation 
(full tank) from the AM condensate vent; hence, the angular 
momentum associated with this vent is negligible. 

4.3 AM Molecular Sieves 


Thrust is calculated by the method of Reference 2 
and is essentially constant over the mission. Thrust imbalance 
for the pair of vents is taken to be 5 percent of the total 
flow,* and impingement on the MDA, OWS and solar arrays is 
calculated in Appendix C. Imbalance and impingement together 
give rise to the numbers shown in Table II. 


♦Obtained from R. E. Tinius, MSFC, October 21, 1970. 
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4 . 4 OWS Urine Pump System 

Four methods of dumping urine into the LOX tank 
are currently under consideration: 

a) dump as an unconfined liquid, 

b) dump in bags that will burst on entering tank, 

c) dump in bags having filters that will pass gas 
but not liquid, and 

d) dump in sealed bags that will not burst. 

Methods a and b are treated here because they represent the 
most severe case. Behavior of the filters in method c is 
poorly understood, and method d of course results in no venting. 

Thirteen pounds of urine are dumped each day. Flow 
rate is obtained from Reference 3, enc. 7, by assuming flow 
rate to be proportional to tank pressure, and specific impulse 
calculations are covered in Appendix D. Three factors influence 
the net propulsive effect of the LOX tank vents, flow imbalance, 
nozzle misalignment, and OWS solar array impingement. Flow 
imbalance, due mainly to the fact that one leg of the plumbing 
is about 30 inches longer than the other, amounts at most to 
about five percent.* Misalignment error calculations are covered 
in Appendix E, and impingement force coefficients are taken from 
Table I of Reference 4. The three effects are summarized below 
in terms of torque per unit thrust (ft- lb/lb) . 




V F 

T /F 

y 

T z /F 

imbalance 


.035 

-2.059 

.917 

misalignment 


+ .045 

+ .377 

+ .826 

impingement 


-1.588 

.225 

- .590 

total (upper 

signs) 

-1.508 

-1.457 

1.153 

total (lower 

signs) 

-1.598 

-2.211 

- .499 


The values for the lower signs are used, since they appear 
to be worst case. 


*This information supplied by Larry Lannon, McDonnell 
Douglas-West, August 10, 1970. 
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Values shown in Table II, which are based on a 
single 13 lb dump, are for the orbit in which the dump 
occurs; venting decreases to zero over the succeeding five 
orbits are shown in the figure. 

4 . 5 OWS Trash Dumps 

Present estimates are that there will be five 
trash dumps per day in perforated bags; one bag will 
contain 4 lb of water, and each of the others will contain 
2 lb of water.* Calculations are similar to those for the 
urine dump, except that flow rate is based on enc. 8 of 
Reference 3. Venting requires two orbits, and values for 
the first orbit of the 4 lb dump are shown in Table II. 

4 . 6 M092 Lower Body Negative Pressure Experiment 

Flow rate is estimated to be 2 cfm,** venting lasts 
15 minutes, and the experiment is assumed to be performed every 
other day. As the vent is located under the OWS meteoroid 
shield, diffusion acts to reduce the net effective thrust. 
Worst-case values in Table II are for flow unimpeded by the 
meteoroid shield, and typical values in the table are simply 
one-third of the worst-case values.*** 

4 . 7 M479 Zero -G Flammability Experiment 

There are six experiment sessions, which we assume 
to occur on different days. Each session comprises twelve 
short, equal-magnitude ventings spaced 5 minutes apart; the 
total accumulation for one session is as shown in Table II. 

5 . 0 Momentum Profiles for Typical and Worst-Case Days 

Here we synthesize composite 24-hour momentum pro- 
files from the individual time histories just discussed. Such 
a synthesis is only an estimate because mission timelines have 
not been established yet. Therefore, these profiles just 
serve as examples of what might occur. 


* Information from R. E. Tinius , MSFC, October 21, 1970. 

** "Experiment Requirements Document for Inflight Lower 
Body Negative Pressure," MSC-KW-D-69-23 REVA, September 1970, 
paragraph 4.1.6b. 

***0ne-third factor suggested by R. E. Tinius, MSFC, 
November 12, 1970. 
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VENTING TORQUE DUE TO A 13 LB 
URINE DUMP (BASED ON REF. 3. ENC. 7) 


Two profiles are presented, one called worst 
case, in which options are chosen and scheduling is done 
so that the total daily momentum accumulation is large 
and much of it is concentrated in one orbit, and one called 
typical, in which decisions are made and scheduling is done 
with a view toward reducing total accumulation and keeping 
the rate of accumulation uniform. No significance beyond 
this should be ascribed to "worst case" or "typical". Table 
III lists specific differences between the two profiles. 

Results for the two profiles are given in Table IV. 




TABLE III 
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TABLE IV 

COMPOSITE ANGULAR MOMENTUM ACCUMULATIONS 
(FT-LB-SEC/ORBIT) 



TYPICAL DAY 

WORST-CASE DAY 1 

ORBIT 

MAG. 

DIR. 

COSINES 

MAG. 

DIR. 

COSINES | 

1 

370 


—.978 

-.205 

7290 

-.145 

-.944 

-.297 

2 

■ fl 


-.167 

.985 

920 

-.547 

-.813 

-.199 

3 



-.167 

.985 

670 

-.544 

-.819 

-.206 

4 

190 

.021 

-.167 

.985 

850 

-.546 

-.815 

-.201 

5 

190 

.021 

-.167 

.985 

360 

-.498 

-.836 

-.229 

6 

190 

.021 

-.167 

.985 

110 

-.300 

-.894 

-.332 

7 

190 

.021 

-.167 

.985 

420 

-.510 

-.831 

-.222 

8 

190 

.021 

-.167 

.985 

130 

-.346 

-.886 

-.310 

9 

190 

.021 

-.167 

.985 

40 


-.749 

-.524 

10 

190 

.021 

-.167 

.985 

420 


-.831 

-.222 

11 

190 

.021 

-.167 

.985 

130 

-.346 

-.886 

-.310 

12 

190 

.021 

-.167 

.985 

40 

3 

-.749 

-.524 

13 

190 

.021 

-.167 

.985 

420 

■ 

-.831 

-.222 

14 

190 

.021 

-.167 

.985 

130 

-.346 

-.886 

-.310 

15* 

190 

.021 

-.167 

.985 

40 


-.749 

-.524 


•THERE ARE ACTUALLY ABOUT 15.4 ORBITS/DAY, BUT THE REMAINDER 
IS NEGLECTED 
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APPENDIX A 


Thrust and Flow Rate for Liquid Venting 


The thrust F is 


F 



A e P e 


(A-l) 


where m = w/g is the mass flow rate, V is velocity, A is 
cross-sectional area, and p is pressure; subscript e 
refers to conditions at the exit plane. From continuity 


m = pAV 
e e 


(A- 2) 


where c denotes mass density. If V c = 0 in the reservoir, 
we get from the Bernoulli equation 


= P< 



(A-3) 


From these relationships 

F = 2A e (p 0 - p e ) + A e p e (A-4) 

= A e (2p 0 - p e ) (A-5) 

Observe that p g lies between p Q and zero in value, probably 
closer to zero. Thus, approximately 
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If P e << Po we see by comparison of (A-l) and (A-4) 

that the first term on the right hand side dominates the 
second term. In this case, approximately 


F - 



(A-7) 


and by use of (A -2) 


• 2 

w 

27 

09 


(A-8) 


With this same assumption that p e <<p 0 , we obtain 
from (A- 2) and (A- 3) that, approximately 
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APPENDIX B 


Minimum Flow Rate for SM Cryogenic Storage Tanks 


HYDROGEN TANKS 


The specified heat leak is 7.25 Btu/hr per tank 
at 140 ° F . The tanks' environment in Sky lab, however, will 
run between -3 and +16°F,* so we need the heat leak at these 
temperatures. If heat transfer is by conduction 


g x - K l( T e - t h > 

and if by radiation 

q 2 = K 2 (T e 4 - T h 4 ) 

where are the heat leak rates, are constants, T g is 
the environmental temperature, and T^ is the temperature 
of the contents. Using the specified heat leak and T = 64° R 
(Reference 5, Figure 8-7) , one can find and K 2 > Then one 
can compute and q 2 as functions of quantity by using 
the T^ vs. quantity plot of Reference 5, Figure 8-7 
(p H = 260 psia) . It turns out that q^ and q 2 are nearlv 
independent of quantity, and q^ = 5.5 Btu/hr, q 2 = 2.8 Btu/hr. 
Since radiation is predominant in this application, we take 
q = 4.0 Btu/hr. 


( B-l ) 


( B-2 ) 


*Obtained from Ray Zajac, North American Rockwell, 
August 24, 1970. 
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Flow Rate vs. Quantity 

We desire the flow rate w = dw/dt given the heat 

leak q = dq/dt. The missing factor, dq/dw, which is plotted 
vs. quantity in Figure B-l, comes from Reference 5, Figure 
8-6a.* Figure B-2 is obtained by dividing q by the values in 
Figure B-l. 



FIGURE B-1 - dq/dw VS. QUANTITY 


*The symbol m for mass is conventional in this work, but 
since lb units are used, we switch to w to maintain consistency 
with other parts of the memorandum. 
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H 2 QUANTITY (W) - LB 


Flow Rate vs. Time 


FIGURE B-2 - FLOW RATE VS. QUANTITY 


This is obtained by integrating the reciprocal of 
Figure B-2 backward in quantity to get a plot of time vs. 
quantity. The boundary condition used is that t = 0 when 
quantity w = 12 lb. This condition corresponds to the nominal 
mission, which with respect to venting is also worst case 
because, as can be seen from Figure B-2, any slippage in the 
schedule will result in a smaller initial value of w and a 

corresponding smaller peak value of w.* 

By use of the time vs. quantity relation just 
described, the abscissa of Figure B-2 can be relabeled to 
give the desired flow rate vs. time plot shown in Figure B-3. 
Note that this figure gives the total flow from both tanks 


♦This information obtained from Walter Scott, MSC, 
October 29, 1970. 



Bellcomm. Inc. 



DAYS AFTER FUEL CELL SHUTDOWN 


FIGURE B 3 - H 2 FLOW RATE VS. TIME 


OXYGEN TANKS 


The calculations for these tanks proceed in the 
same way as those for the hydrogen tanks just discussed, 
and the figures are based on an intital quantity of 112.5 
lb per tank for the nominal mission.* Loss of 0 2 from the 


cabin atmosphere due to leakage and metabolic needs is 
estimated at 0.55 lb /hr.* Figure B-4, from which venting 
thrust is calculated, shows the flow rate from both tanks 
in excess of the 0.55 lb/hr. This excess may either be 
vented overboard or it may be added to the cabin atmosphere 
in which case it will raise the 0 2 partial pressure by 0.24 


9 


psi. If there is less than 62 lb left in each tank at fuel 
cell shutdown, there will be no excess flow. 


♦Obtained from Mr. Scott. 
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APPENDIX C 


Molecular Sieve Impingement Analysis 


The flow field is calculated by the method of 
Hill and Draper (Reference 6), which has been suggested 
by NASA (see Reference 7). Specifically; Hill and Draper 
propose the following formula^for distribution of the mass 
flux per unit solid angle, d m/dft, as a function of 0, 
the angle off of the vent axis. 


(dm/dft) 0 
(dm/dfi) e=Q 


(C-l) 


Sufficiently far from the vent, X approaches the value 


‘ = [ V” (1-C /C )] 

max 


(C-2) 


The thrust coefficient ratio C^/C^ can be obtained from 

F F max 

equations (10) and (30) of Reference 8. For a mean specific 
heat ratio v=1.34 (Reference 1, p. 71) and unit area ratio 


C^/C^ = 0.6656 

A max 

and 

X = 1.69 

The unknown factor (dm/dfi) Q _ 0 in (C-l) can by 
determined by using (C-l) in a continuity relationship 
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M = 


/ 


dm 

do 


dft 


= 2 it (dm/dO) 


6=0 




(1-COS0) 


2 


sine 


de 


The integral in 6 above can be reduced to the error integral, 
which is tabulated in Reference 8. The result is 


(dm/dn) 0=o 



(C-3 ) 


Since impingement is nearly symmetrical in the +y 
and -y directions, the impingement torque is essentially 
parallel to the y axis. Thus, we compute just the torque 
T=RF about the y axis, where R is the effective lever arm 
and F is the impingement force on a surface of solid angle 
For perfectly elastic collisions 


F = 2e ^ VV (C-4) 


where e is the cosine of the angle between the streamline 
and the surface normal and where V is the flow speed. Use 
of (C-l) , (C-3) and (C-4) in T=RF yields 


T 


KeRfie 


•x 2 (l-cose) 2 


(C-5) 


where 

K = - MV A 

7T 

-3 2 

has the value 4.56 x 10 slug-ft/sec for the molecular 
sieve vents. The approach is now straightforward: 



BELLCOMM, Inc. 


C-3 


determine which surfaces the molecular sieve effluent 
can impinge upon; 

divide these surfaces into sufficiently small subsur- 
faces ; 

calculate t, R, fi, and 0 for each subsurface; 

calculate T for each subsurface using (C-5) ; and 

add the individual torques to find the total impingement 
torque . 

This approach was used on the following surfaces: 

OWS solar arrays (mostly shadowed by the OWS) ; 

ATM solar arrays (mostly shadowed by the OWS, STS, 
and MDA) ; 

STS and MDA (broken into 24 subsurfaces) ; and 

AM truss meteoroid curtain (broken into 10 subsurfaces) . 

In each case impingement torques were computed for only one 
vent of the pair, and the results were doubled (this is 
accounted for in the numerical value of K) . The resulting 

_ 3 

total torque is 6 x 10 ft- lb, the major portion of which 

comes from the STS/MDA surfaces. 
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APPENDIX D 

Specific Impulse for the LOX Tank Vents 


By use of a specific heat ratio y=1-28 for water 
and an area ratio A e /A t = (4/1. 215) 2 for the nozzle, one 
obtains from Reference 8 (eqn (19) , p. 67) a pressure ratio 



.009015 


where e denotes conditions at exit and c denotes conditions 
in the tank (chamber) .* Use of this pressure ratio and the 
value T c = 510°R** in formula (16), p. 62, of Reference 8 

gives 


V = 2878 ft/sec 
e 


Furthermore, one obtains from eqn (13), p. 60, of the same 
reference 


P c = 2.23 x 10 5 m lb/ft 2 

where m is the nozzle mass flow rate. It follows from the 

A 

above pressure ratio that p g =201 m lb/ft , and from eqn (A-l) 


* Mean free path calculations for this vent indicate that 
continuum theory is valid here. 

**Suggested by D. P. Woodard. 
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D-2 


I 


_F 

mg 



A p 
e^e 

mg 


89. 


sec 


I is the specific impulse of one 


LOX tank nozzle 
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APPENDIX E 


Misalignment of "Nonpropulsive" Vents 


Nominally, the OWS "nonpropulsive" vents are mounted 
on opposite sides of the cylindrical Workshop so that their 
lines of action coincide. The objective here is to determine 
what torque is apt to result if the mounting varies from 
nominal within specified tolerances. 

As shown in the figure, the first vent of the pair 
is nominally mounted at a distance R from the Workshop longi- 
tudinal (x) axis and at an angle 6 about this axis, and the 
thrust vector is parallel to the location vector. Let sub- 
scripted Greek letters denote small random deviations from 
nominal. 6^ is the angular mounting deviation, is the 

axial mounting deviation, and and are the angular 

deviations of the thrust vector relative to the location 
vector, these being right-hand rotations about x and z, 
respectively, when 6 = 0 . Mounting of the second vent is 
described similarly, except that the angular location 6+62 

is measured from -y toward -z. 


z 



The location vectors relative to the mass center 

(x c' V z c> are 
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r. = (e.-x ) i+ t+Rcos (0 + 6 . ) -y ] j_+ [+Rsin ( 6 + 6 . ) -z ] k 

X 1 C — 1 C ~ 1 C ” 


(E-l) 


the upper signs referring to the first vent (i=l) and the 
lower ones referring to the second vent (i=2) . The thrust 
vectors are 


F^=+F [ c ^i^-cos (e+6 i +C i )j_-sin (e + 6 i + C i >k] 


(E-2) 


and the torque is 


T = £]_ X £]_ + £2 X — 2 


( E— 3 ) 


If (E-l) and (E-2) are substituted into (E-3) and 
quadratic and higher terms in the small deviations are dropped, 
we obtain 


T=F{ [ (y c C+z c S) (6 1 -6 2 )-(R-y c C-Z c S)C 1 -(R+y c C+Z c S) C 2 U 
+ [-x c S ( 6 6 2 ) +S ( £ ]_“ e 2 ) _x c c ( C i“ C 2 > "*" ( RS ~ Z C ) C i+ (RS+z c ) 

+ [-x c S(6 1 -6 2 )-C(E 1 -E 2 )-x c S(e 1 -C 2 )-(RC-y c ) ? 1 ~(RC-y c ) ^]k} (E-4) 

where C = cose and S = sine. 

Let T , T , T be the i , j , k measure numbers of 
x y z 

T, and let 



T 


(E-5) 
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We desire the covariance, p, of T, 


P = E{ (T-T q ) (T-T Q ) ' } 


(E-6) 


where E is the expectation operator, where Tq is the mean of 

T, and where primes denote transposition. With the assumption 
that all of the deviations are zero mean and uncorrelated, the 
elements of P may be computed by inspection from (E-4) . For 
example , 


p ll - E!T x 2! ' F2 i ( yc c+z c s > 2(o «i 2+< ’ 5 2 2 > 

+ ' R -y c c - z c s) 2 °u 2+ (R+ y c c+ z o s) °c2 2 ' (E - 7) 

2 2 

where =E{6^ } is the variance of 6^, etc. If the 

deviations are uniformly distributed, as we shall assume, 

2 o 

z = A /3 , where +A is the tolerance on the deviation under 
consideration. 

Application to LOX Tank Vents 

Values for x , y , z , R, and 0 come from Tables I 
c c c 

and II. The tolerances are* 


e 1 , z 2 : + 0.5 inch 

6 x , § 2 , iy K 2 ’ ^2 : ± 0 * 5 de 9 ree 

For these values 


*Supplied by George Carlson, McDonnell Douglas-West, 
August 10, 1970. 
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1 




P 


0064 

.0057 

.0083 

0057 

.0416 

.0751 

0083 

.0751 

.1719 


With a view toward finding a coordinate transformation that 
will eliminate the substantial off-diagonal elements of P, 
eigenvalues and eigenvectors of P are calculated. The eigen- 
values are 


. 0047F 2 , .0085F 2 , 0.2067F 2 


Comparison of the relative magnitudes of these three numbers 
reveals that the torque uncertainty in the directions of the 
first two eigenvectors is much less than it is in the direction 
of the third eigenvector. Thus, in all likelihood, the mis- 
alignment torque will be closely aligned with the third eigen- 
vector, and this gives the direction of the misalignment torque. 
The magnitude is chosen so that there will be about 95% prob- 
ability that the actual torque will be smaller than the value 
calculated. This is done by invoking the central limit 
theorem, which says that misalignment torque will tend toward 
a Gaussian distribution even if the deviations are not 
Gaussian, by taking the square root of the third eigenvalue, 
which gives standard deviation a, and by multiplying by two, 
since 2o gives about 95% assurance that the actual value will 
be less than the calculated value. The resulting values for 
misalignment torque are 


+.045F, 


+.377F, 


+. 826F 
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